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ABSTRACT 



We study a general quantum system interacting with environment 
modeled by the bosonic heat bath of Caldeira and Leggett type. General 
interaction Hamiltonians are considered that commute with the system's 
Hamiltonian so that there is no energy exchange between the system and 
bath. We argue that this model provides an appropriate description of 
adiabatic quantum decoherence, i.e., loss of entanglement on time scales 
short compared to those of thermal relaxation processes associated with 
energy exchange with the bath. The interaction Hamiltonian is then 
proportional to a conserved "pointer observable." Calculation of the ele- 
ments of the reduced density matrix of the system is carried out exactly, 
and time-dependence of decoherence is identified, similar to recent re- 
sults for related models. Our key finding is that the decoherence process 
is controlled by spectral properties of the interaction rather than system's 
Hamiltonian. 
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1. Introduction 



Quantum decoherence, dissipation, and thermalization due to inter- 
actions with environment have long been important fundamental issues 
theoretically and experimentally. Decoherence and related topics 

have attracted much interest recently due to rapid development of new 
fields such as quantum computing and quantum information theory/^^"^^) 
Decoherence due to external interactions is a major obstacle in the way of 
implementation of devices such as quantum computers. Thus in addition 
to studies of the physics of decoherence processes there emerged a new 
field of quantum error correction^^^"^^) aiming at effective stabilization 
of quantum states against decoherence essentially by involving many ad- 
ditional quantum systems and utilizing redundancy. The present work 
contributes to the former topic: the physics of decoherence. 

Decoherence is a result of the coupling of the quantum system under 
consideration to the environment which, generally, is the rest of the uni- 
verse. In various experimentally relevant situations the interaction of the 
quantum system with environment is dominated by the system's micro- 
scopic surroundings. For example, the dominant source of such interac- 
tion for an atom in an electromagnetic cavity is the electromagnetic field 
itself coupled to the dipole moment of the atom.*^^^^ In case of Josephson 
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junction in a magnetic flux*^^^) or defect propagation in solids, the interac- 
tion can be dominated by acoustic phonons or delocalized electrons. 
Magnetic macromolecules interact with the surrounding spin environ- 
ment such as nuclear spins. Numerous other specific examples could 
be cited. 

In this work we aim at a general phenomenological description that 
models the physically important effects of external interactions as far as 
adiabatic decoherence, to be defined later, is concerned. We note that 
generally thermalization and decoherence are associated with the inter- 
action of the quantum system, described in isolation by the Hamiltonian 
Hs, with another, large system which we will term the "bath" and which 
internally has the Hamiltonian Hb- The actual interaction will be rep- 
resented by the Hamiltonian Hj so that the total Hamiltonian of the 
system, H, is 

H = Hs + Hb + Hi. (1.1) 

It is important to realize that typically the bath is a large, macroscopic 
system. Truly irreversible interactions of a quantum system with its en- 
vironment, such as thermal equilibration or decoherence associated with 
measurement processes, can only be obtained in the Hamiltonian descrip- 
tion (1.1) when it is supplemented by taking the limit of the number of 
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particles or degrees or freedom of the bath going to infinity. 

Interactions of a quantum system with macroscopic systems can lead 
to different outcomes. For instance, interaction with a true "heat bath" 
leads to thermalization: the reduced density matrix of the system ap- 
proaches exp {—f3Hs) for large times. Here 

/? = 1/ikT) (1.2) 

as usual, and by "reduced" we mean the density matrix traced over 
the states of the bath. On the other hand for decoherence we expect 
the reduced density matrix to approach a diagonal form in the "pre- 
ferred basis" somehow selected by the "pointer observable" Hermitian 
operator which is thereby "measured" by the macroscopic sys- 
tem (bath). 

It is important to realize that study of decoherence in the present 
context does not fully resolve the problem of understanding quantum 
measurement and other fundamental issues at the borderline of quantum 
and classical behaviors, such as, for instance, the absence of macroscopic 
manifestations of Schrodinger-cat type quantum superposition of states. 
The more optimistic recent literature^^^^-* considers description of entan- 
glement and decoherence the key to such understanding. However, these 
fundamental problems have remained open thus far. 
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The most explored and probably most tractable approach to model- 
ing the environmental interactions has involved representing their effects 
by coupling the original quantum system to a set of noninteracting har- 
monic oscillators (bosonic heat bath).'^^'^'^~^^'^^'^^~^^^ Fermionic heat 
bath can be also considered, e.g., Ref. 34. We will use the term "heat 
bath" for such systems even when they are used for other than thermal- 
ization studies because they have the temperature parameter defined via 
initial conditions, as described later. 

Rigorous formulation of the bosonic heat bath approach was initi- 
ated by Ford, Kac and Mazur^^^^ and more recently by Caldeira and 
Leggett.'^^^'^^-' It has been established, for harmonic quantum systems, 
that the influence of the heat bath described by the oscillators is effec- 
tively identical to the external uncorrelated random force acting on a 
quantum system under consideration. In order for the system to satisfy 
equation of motion with a linear dissipation term in the classical limit 
the coupling was chosen to be linear in coordinates while the coupling 
constants entered lumped in a spectral function which was assumed to 
be of a power-law form in the oscillator frequency, with the appropriate 
Debye cutoff. We will make this concept more explicit later. 

This model of a heat bath was applied to studying effects of dis- 
sipation on the probability of quantum tunneling from a metastable 
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state. "^^'^^^ It was found that coupling a quantum system to the heat bath 
actually decreases the quantum tunneling rate. The problem of a parti- 
cle in a double well potential was also considered. In this case the 
interaction with the bath leads to quantum coherence loss and complete 
localization at zero temperature. This study has lead to the spin-boson 
Hamiltonian^^'^^) which found numerous other applications. The Hilbert 
space of the quantum systems studied was effectively restricted to the 
two-dimensional space corresponding to the two lowest energy levels. 

Another possible application of the bosonic heat bath model con- 
cerns aspects of quantum measurement. It is believed that the bath is 
an intrinsic part of a measuring device. In other words, it continuously 
monitors the physical quantity whose operator is coupled to it.^^"^) It 
has been shown in the exactly solvable model of the quantum oscillator 
coupled to a heat bath*^^^ that the reduced density matrix of the quan- 
tum system decoheres, i.e., looses its off-diagonal elements representing 
the quantum correlations in the system, in the eigenbasis of the interac- 
tion Hamiltonian. It has also been argued that the time scale on which 
this "measurement" occurs is much less than the characteristic time for 
thermal relaxation of the system. 

It is natural to assume that if such a "bath" description of the pro- 
cess of measurement of a Hermitian operator As exists, then the inter- 
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action Hamiltonian Hj in (1.1) will involve A5 as well as some bath- 
Hilbert-space operators. No general description of this process exists. 
Furthermore, when we are limited to specific models in order to obtain 
tractable, e.g., analytically solvable, examples, then there is no general 
way to separate decoherence and thermalization effects. We note that 
thermalization is naturally associated with exchange of energy between 
the quantum system and heat bath. Model system results and general 
expectations mentioned earlier suggest that at least in some cases de- 
coherence involves its own time scales which are shorter than those of 
approach to thermal equilibrium. 

In this work we propose to study adiabatic decoherence , i.e., a special 
case of no energy exchange between the system and bath. Thus we assume 
that Hs is conserved, i.e., [Hs^H] = 0. This assumption is a special 
case of "quantum nondemolition measurement" concept '-^'^^^ exemplified 
by the Kerr effect, for instance. Since Hs and Hb is (1.1) operate in 
different Hilbert spaces, this is equivalent to requiring 



[Hs,Hi]=0. 



(1.3) 



Furthermore, we will assume that Hj is linear in As: 



Hi = AsPb 



(1.4) 
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where Pb acts in the Hilbert space of the bath. Then we have 

[As,Hs] = 0. (1.5) 

Thus, we consider cases in which the measured, "pointer" observable As is 
one of the conserved quantities of the quantum system when it is isolated. 
Interaction with the bath will then correspond to measurement of such 
an observable, which can be the energy itself. Specifically, the model 
of Ref. 14 corresponds to A^ = Hs for the case of the spin-| two-state 
system, motivated by quantum-computing applications; see also Refs. 2, 
12-15. The models of Refs. 1 and 2 correspond to the choices of A^' = Hs 
and As- = f{Hs), respectively, for a system coupled to a bosonic spin 
bath, where / is an arbitrary well-behaved function. 

Here we derive exact results for adiabatic decoherence due to cou- 
pling to the bosonic heat bath, assuming general As that commutes with 
Hs- While technically this represents an extension of the results of Refs. 1 
and 2, we demonstrate that the general case reveals certain new aspects of 
the decoherence process. Our new exact-solution method utilizes coher- 
ent states and may be of interest in other applications as well. In Section 
2, we define the system. Specifically, we choose the bosonic heat bath 
form for Hb and Pb in (1.1) and (1.4), but we keep Hs and A^ general. 
However, we also analyze the mechanism leading to exact solvability of 
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general models of this type. Section 3 reports our derivation of the exact 
expression for the reduced density matrix of the system. Discussion of 
the results and definition of the continuum limit are given in Section 4. 
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2. Models of Adiabatic Decoherence 

We will be mainly interested in the following Hamiltonian for the 
quantum system coupled to a bath of bosons (harmonic oscillators) la- 
beled by the subscript k: 

H = Hs + uJkoiak + {glo^k + QkOk) ■ (^-l) 

k k 

f 

Here and are bosonic creation and annihilation operators, respec- 
tively, so that their commutation relation is [a^, a|] = 1. The second term 
in (2.1) represents the free field or Hamiltonian of the heat bath Hb- The 
last term is the interaction Hamiltonian Hj. The coupling constants will 
be specified later; exact results obtained in Section 3 apply for general 
Uk and Qk- Here and in the following we use the convention 

n = 1 (2.2) 

and we also assume that the energy levels of each oscillator are shifted 
by ^ujk so that the ground state of each oscillator has zero energy. 

Since we assume that Hs and commute, we can select a common 
set \i) of eigenstates: 



Hs\i) = Ei\i) 



(2.3) 



As|i) = Aj|i> . 



(2.4) 



One of the simplifications here, due to the fact that Hs and commute, 
is that these eigenstates automatically constitute the "preferred basis" 
mentioned earlier. 

We will assume that initially the quantum system is in a pure or 
mixed state described by the density matrix p{0), not entangled with the 
bath. For the bath, we assume that each oscillator is independently ther- 
malized (possibly by prior contact with a "true" heat bath) at tempera- 
ture T, with the density matrix 9k. The total system-plus-bath density 
matrix will then be the product 



(2.5) 



k 



Here 




(2.6) 
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Z, = (l-e-''"')-i, 



(2.7) 



where is the partition function for the oscillator k. The quantity j3 
was defined in (1.2). Introduction of the temperature parameter via the 
initial state of the bath is common in the literature/^'^'^"^^'^^"^^'^^'^^'^^^ 
While it may seem artificial, we recall that the bath is supposed to be 
a large system presumably remaining thermalized on the time scales of 
interest. Specific results indicating that the bosonic heat bath can be 
viewed as a source of thermalizing noise have been mentioned earlier; see 
also Ref. 35. 

Our objective is to study the reduced density matrix of the system 
at time t > 0; it has the following matrix elements in the preferred basis: 



Here the outer trace is taken over the states of the heat bath, i.e., the 
bosonic modes. The inner matrix element is in the space of the quantum 
system. Note that for no coupling to the bath, i.e., for — 0, the density 
matrix of the system is simply 




-iHt 






gk=o 



(2.9) 
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For the interacting system, the heat-bath states must be summed 
over in the trace in (2.8). It is instructive to consider a more general case 
with the bath consisting of independent "modes" with the Hamiltonians 
Mfe, so that 

HB = J2Mk, (2.10) 

k 

f 

where for the bosonic bath we have = ujk%ak- Similarly, for the 
interaction term we assume coupling to each mode independently, 

Hi = AsY,Jk, (2.11) 

k 

where for the bosonic bath we have Jk = g^^k + 90^- Relation (2.5) 
remains unchanged, with the definitions (2.6) and (2.7) replaced by 

Ou = Z-^e-''^'^ , (2.12) 

Zk = Trfc [e-^^'^] , (2.13) 

where the trace is over a single mode k. 

Owing to the fact that Hs and As share common eigenfunctions, 
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the inner matrix element calculation in (2.8), in the system space, can be 
expressed in terms of the eigenvalues defined in (2.3)-(2.4). Specifically, 
we define the bath-space operators 



(2.14) 



k k 

which follow from the form of the Hamiltonian. The calculation in (2.8) 
then reduces to 



Pmn{t) = Ttb 



m\e 



k 



which yields the expression 



(2.15) 



Pmn{t) = Pmn (O)TrB 



k e 



ih„ t 



(2.16) 



We will now assume that the operators of different modes k com- 
mute. This is obvious for the bosonic or spin baths and must be checked 
explicitly if one uses the present formulation for a fermionic bath. Then 
we can factor the expression (2.16) as follows: 



(2.17) 



This expression, or variants derived in earlier works, ^'^'^'^ suggests that 
the problem is exactly solvable in some cases. Indeed, the inner trace in 
over a single mode of the bath. For a spin bath of spin-| "modes" the 
calculation involves only (2 x 2)-matrix manipulations and is therefore 
straightforward. '^'^'^ However, in this case the only nontrivial choice of 
the "pointer observable" corresponds, in our notation, to As = Hs, with 
both operators usually chosen equal to the Pauli matrix cr^ . There is also 
hope for obtaining analytical results for other baths with modes in finite- 
dimensional spaces, such as spins other than ^ ; we have not explored this 
possibility. 

For the bosonic spin bath, the calculation is in the space of a single 
harmonic oscillator. It can be carried out by using operator identities. ^'^ 
We have used instead a method based on the coherent-state formalism 
which is detailed in the next section. 
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3. Exact Solution for the Density Matrix 



We utilize the coherent-state formaUsm, e.g., Refs. 35, 36. The co- 
herent states \z) are the eigenstates of the annihilation operator a with 
complex eigenvalues z. Note that from now on we omit the oscillator 
index k whenever this leads to no confusion. These states are not orthog- 
onal: 



They form an over-complete set, and one can show that the identity 
operator in a single-oscillator space can be obtained as the integral 



{zi\z2) = exp zlz2 - - \zi 




) 



(3.1) 




(3.2) 



Here the integration by definition corresponds to 



(f'z = —d (Rez) d (Imz) . 



(3.3) 



TT 



Furthermore, for an arbitrary operator we have, in a single-oscillator 

space. 



TtA = J (fz{z\A\z). (3.4) 
Finally, we note the following identity/^^^ which will be used later, 

In this relation is an arbitrary c-number, while M denotes normal 
ordering. 

The result (2.17) for the reduced density matrix, assuming the bosonic 
spin bath, can be written as 



(3.5) 



Pmnif) — Pmn(0)e^ " '^^ J_ J_ Smn,k — [Pmn(^)]g^=o J. J. ^rnn,k 5 (3-6) 

k k 

where we used (2.9). Omitting the mode index k for simplicity, the 
expression for Smn for each mode in the product is 



Smn — Z Tr 



(3.7) 



where the trace is in the space of that mode, and we defined 



7m = ujo^a + Am (g*a + ga^^ . 



(3.8) 



The partition function Z is given in (2.7). Relations (3.6)-(3.8) already 
illustrate one of our main results: apart from the phase factor which 
would be present in the noninteracting case anyway, the system energy 
eigenvalues do not enter in the expression for Pmn{t)- The interest- 
ing time dependence is controlled by the eigenvalues of the "pointer 
observable" operator As (and by the heat-bath coupling parameters ojk 
and Qk). 

In order to evaluate the trace in (3.7), we use the coherent-state 
approach. We have 



ZS„„ = J Szf,SziSz2 (2o|e""^"'ki)(2i|e"''"''^°k2)(22|e"^"|2o>- 

(3.9) 

The normal-ordering formula (3.5) then yields for the middle term. 



exp 



(3.10) 



In order to evaluate the first and last factors in (3.9) we define shifted 
operators 
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(3.11) 



in terms of which we have 

Jm = u;r,U-Xl^uJ-'\9f■ (3.12) 

Since 77 and 77^ still satisfy the bosonic commutation relation [77,77^] = 1, 
the normal-ordering formula applies. Thus, for the first factor in (3.9), 
for instance, we get 

U|e-''^-ki> = e«^™'#Uki)e('="'"'-')(^°-'''"'^)<^'^''"*' . (3.13) 

Collecting all these expressions, one concludes that the calculation of 
Smn involves six Gaussian integrations over the real and imaginary parts 
of the variables zo,zi,Z2. This is a rather lengthy calculation but it can 
be carried out in closed form. The result, with indices k restored, is 

Smn,k = exp {-UJ^^\gk\'^Pmn,k) , (3.14) 

where 
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Pmn,k = 2 (A^ - Xnf sin^ ^ coth + i (A^ - A^) (sina;^^ - Ukt) . 

(3.15) 

The expression (3.15), with (3.14), when inserted in (3.6), is the prin- 
cipal result of this section. It will be discussed in the next section. Here 
we note that in the studies of systems involving the bosonic heat bath 
one frequently adds the "renormalization" term^'^^ in the Hamiltonian, 



H = Hs + Hb + Hi + Hr, (3.16) 

where in our case 



HR = AlJ2u;^'\gu\\ (3.17) 

k 

The role of this renormalization has been reviewed in Ref. 29. Here we 
only notice that the sole effect of adding this term in our calculation is to 
modify the imaginary part of Pmn,k which plays no role in our subsequent 
discussion. The modified expression is 



^mn,k = 2 (A^ - Xnf siu^ coth + i (A^ - A^) sin ujkt . (3.18) 



4. Continuum Limit and Discussion 



The results of the preceding section, (3.6), (3.14), (3.15), can be 
conveniently discussed if we consider magnitudes of the matrix elements 
of the reduced density matrix p{t). We have 



\Pmn{t)\ = |pmn(0)|exp 



liXm-Xnfr{t) 



(4.1) 



where we introduced the factor ^ to have the expression identical to that 
obtained in Ref. 14: 



r(t) =8^c.,-^b|^sin^^coth^. (4.2) 

k 

These results suggest several interesting conclusions. First, the decoher- 
ence is clearly controlled by the interaction with the heat bath rather 
than by the system's Hamiltonian. The eigenvalues of the "pointer ob- 
servable" As determine the rate of decoherence, while the type of the 
bath and coupling controls the form of the function T{t). It is interesting 
to note that states with equal eigenvalues will remain entangled even 
if their energies Em are different. As expected, the magnitude of the 
diagonal matrix elements remains unchanged. 

Secondly, we note that T{t) is a sum of positive terms. However, for 
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true decoherence, i.e., in order for this sum to diverge for large times, one 
needs a continuum of frequencies and interactions with the bath modes 
that are strong enough at low frequencies; see below. From this point on, 
our discussion of the function T{t) is basically identical to that in Ref. 14 
(see also Ref. 1); we only outline the main points. In the continuum limit, 
exemplified for instance by phonon modes in solid state, we introduce 
the density of states G{uj) and sum over frequencies rather than modes 
characterized by their wave vectors. The latter change of the integration 
variable introduces the factor which we will loosely write as it must 
be calculated from the dispersion relation of the bosonic modes. Thus we 
have 



r{t) (X fdu ^G(a;)|^(a;)|2 u'^ sin^ ^ coth ^ . (4.3) 
J dhj 2 2 

In Ref. 14, the following choice was considered, motivated by properties 
of the phonon field in solids; see also Refs. 8-11, 12-18, 29: 



^G(a;)|p(a;)p oc uj^'e'^c . (4.4) 



This combination of the coupling constants and frequencies has been 
termed the spectral function. Here Uc is the Debye cutoff frequency. 

Specifically, the authors of Ref. 14 have analyzed the cases n = 1 and 
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n = 3. For n = 1, three regimes were identified, defined by the time scale 
for thermal decoherence, /3, which is large for low temperatures, see (1.2), 
and the time scale for quantum-fluctuation effects, o;^^. Recall that we 
use the units h = 1. The present treatment only makes sense provided 
uj~^ <^ (3. According to Ref. 14, the first, "quiet" regime t <^ uj~^ 
corresponds to no significant decoherence and V oc {uctY. The next, 
"quantum" regime, o;"^ <^ t <^ ^, corresponds to decoherence driven 
by quantum fiuctuations and F oc In(cjc^)- Finally, for t ^ /3, in the 
"thermal" regime, thermal fiuctuations play major role in decoherence 
and T (xt/(3. 

For n = 3, decoherence is incomplete. Indeed, while n must be 
positive for the integral in (4.3) to converge, only for n < 2 we have 
divergent T(t) growing according to a power law for large times (in fact, 
oc t'^~'^) in the "thermal" regime. Thus, strong enough coupling \g{uj)\ to 
the low-frequency modes of the heat bath is crucial for full decoherence. 

In summary, we derived exact results for the model of decoherence 
due to energy-conserving interactions with the bosonic heat bath. We find 
that the spectrum of the "pointer observable" that enters the interaction 
with the bath controls the rate of decoherence. The precise functional 
form of the time dependence is determined both by the choice of heat-bath 
and system-bath coupling. However, for the case studied, it is universal 
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for all pointer observables and for all the matrix elements of the reduced 
density matrix. 
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